Time delay correlations and resonances in ID disordered systems 
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The frequency dependent time delay correlation function K(yi) is studied analytically for a 
particle reflected from a finite one-dimensional disordered system. In the long sample limit K(Q) 
can be used to extract the resonance width distribution p(T). Both quantities are found to decay 
algebraically as F~" , and Q.~ v , v ~ 1.25 in a large range of arguments. Numerical calculations for 
the resonance width distribution in ID non-Hermitian tight-binding model agree reasonably with 
the analytical formulas. 



For a long time one dimensional disordered conduc- 
tors set an example for understanding the localization 
phenomenon in real electronic systems. In this commu- 
nication we consider the reflection of a spinless particle 
from the disordered region of length L and address such 
quantities as the resonance widths distribution and the 
time delay correlation function. 

The issue of time delays and resonances attracted con- 
siderable attention in the domain of chaotic scattering 
jjj and mesoscopics because of relevance for properties 
of small metallic granulas and quantum dots An 
essential progress in understanding of statistical fluctu- 
ations of these and related quantities turned out to be 
possible recently in the context of random matrix 

description of chaotic motion of quantum particles in the 
scattering region. Such a description assumes from the 
very beginning that the time of relaxation on the energy 
shell due to scattering on impurities (Thouless time) is 
negligible in comparison with the typical inverse spac- 
ing between neighboring energy levels (Heisenberg time). 
This assumption neglects effects of Anderson localization 
completely. At the same time the localized states are ex- 
pected to effect the form of the resonance width statistics 
and time delay correlations considerably. 

Difficulties in dealing with the localization analytically 
in higher dimensions suggest one-dimensional samples 
with the white noise random potential as a natural ob- 
ject of study. Various aspects of particle scattering in 
such systems were addressed by different groups recently 
|j|-[?J . Still, the statistical properties of resonances in such 
systems were not studied to our best knowledge. 

It is well-known that the level spacing statistics be- 
comes Poissonian with onset of the Anderson localiza- 
tion when the system size L far exceeds the localiza- 
tion length £ and eigenstates do not overlap any longer. 
Notwithstanding the fact that the argumentation is ap- 
plied to a closed system the same should be true for its 
open counterpart as well. To exploit this fact we de- 
fine the two-dimensional spectral density p(Z) associated 
with broadened levels (resonances) p(Z) = ^2 n S(X — 
ReZ n )S(Y + ImZ n ), where Z = X + iY is the complex 
energy, Z n — E n — iY n /2 is the coresponding position of 
the n-th resonance whose width is r„. We further define 



the resonance correlation function as follows: 

{~ P {z 1 )p{z 2 )) c = (p{z 1 ))s<- 2 \z 1 - z 2 )-y 2 (z 1 ,z 2 ) (i) 

with the brackets standing for the averaging over dis- 
order. The so-called cluster function y 2 {Z\,Z%) reflects 
eigenvalue correlations and therefore should vanish in the 
thermodynamic limit L — > oo. This property provides us 
with a possibility to relate the averaged density of res- 
onances in the complex plane (p(Z)) to the time delay 
correlation function. 

We recall that the time delay at a fixed real energy E 
can be written in terms of p(X, Y) as (see e.g. |H,||) 



t{E) 



2Yp(X, Y) 



Y 2 



(2) 



As far as we restrict ourselves to the thermodynamic limit 
the correlation of time delays at two different energies 
E + £l, and E—£l can be expressed in terms of the corre- 
lation function Eq. ([!]) with the cluster function neglected. 
The formula obtained can be further simplified employ- 
ing the integration over X and taking into account that 
semiclassically Q, Y -C E. In addition we find it to be 
convenient to rescale the spectral density and introduce 
the dimensionless resonance widths distribution function 
p{y) by means of the relation 



(p(X,Y)) = 2™f (X)p(27r^(X) Y) 



(3) 



where v%{E) ~ £,(2n\/E)~ 1 is the density of states cor- 
responding to the system of the size £. This is done, the 
relation between the time delay correlation function and 
the widths distribution takes the form: 



A|(t(£ + uA ( )t(E - wA 5 )) c 

JcT dyz&FPiv) 



(4) 



where = (27r^) _1 standing for the mean level spacing 
corresponding to a single localization volume is a natural 
energy scale for measuring the resonance width as well as 
the time delay correlations. As we will see below Eq.(^) 
gives us a possibility to restore the density p(y) from the 
time delay correlation function. 
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Before we proceed with our analysis it is instructive 
to develop a simple intuitive picture about the form of 
the resonance width distribution for a particle reflection 
from one dimensional disordered sample. We may say 
qualitatively that this distribution (normalized against 
the number of states N) is given by 



iV£ N f L ( 

V{y) = —p{y) + j- J dx 6[y-e~ 



■x/e. 



L 



p(y) + -x (e" i/c < y < e" 1 ) 



(5) 



(6) 



where y = T/(T), and , and x(v) is a characteristic func- 
tion of the interval. The variable x in Eq.(|^) stands for 
the distance from the open edge, and we assume the op- 
posite edge of the sample to be closed for the sake of 
simplicity. The function p(y) is the resonance width dis- 
tribution corresponding to resonances residing at a dis- 
tance x < £ from the open edge. Other states which 
are localized in the bulk give rise to exponentially small 
probability for the particle to escape from the sample 
and the corresponding width is T ~ exp (— This is 
oversimplified picture can give us an idea of y^ 1 depen- 
dence of the resonance density in a rather wide paramet- 
ric range. However our explicit calculations presented 
below demonstrate a slightly different power law: y . 

In the present communication we give only a sketch 
of the calculations, relegating details to a more extended 
publication. To start with, we consider an elastic scatter- 
ing of a spinless particle ip(x) on ID disordered sample. 
The boundary condition ip(0) = is imposed at the origin 
(left edge) and the particle is subject to the white-noise 
potential V(x), (V(x)) — 0, inside the interval (0, L). No 
potential barrier is assumed at the right edge which cor- 
responds to the perfect coupling between the system and 
the continuum. The Shrbdinger equation on the inter- 
val (0, L) can be formulated in terms of the logarithmic 
derivative of the wave function z = dln^/dx plT 



^ = V(x)-E-z 2 
dx 



(7) 



where (V(x)V(y)) — a g 5(x — y), E = k 2 . In what fol- 
lows we scale the variance in such a way that a g = 4E'/£, 
with £ being the localization length and use the rescaled 
potential v(x) = V(x)/2k with the variance l/£. 

After the substitution z = k cot (f> Eq. ([5]) is reduced to 
the equation on the phase <p. The phase shift between 
in- and outgoing waves can be checked to be equal to 2<j>. 
Its derivative over the energy is nothing but the corre- 
sponding time delay. 



k dk 



(8) 



To address time delay correlations at two different ener- 
gies we have to take into account the set of the following 
four equations (cf. WM)- 



■ 1 

4>2 

h 

T-2 



k + q — 2v(x) sin 2 (f>i 
k — q — 2v(x) sin 2 (f>2 
k^ 1 — 2t 1 v(x) sin 2</>i 
fc _1 — 2t2v{x) sin 202 



(9) 



where the dot states for the derivative over x. The first 
two equations are those for phases each related to its own 
eigenstate E = (k + q) 2 , and E — (k — q) 2 . Two more 
equations are those for the time delays. They emerge 
after taking the derivative of the corresponding phases 
over the energy E = k 2 and exploiting an approximation 
q -C k. Eqs. @ contain a multiplicative random force 
vjx) and must be understood in the Stratonovich sense 
§• 

As usual, significant simplifications occur in semiclas- 
sical limit when we restrict ourselves to the large ener- 
gies: k ^> q, kL 3> 1 flip! . Last condition means that 
the joint probability density P(4>i, <f>2,Ti, t%) satisfying a 
Fokker-Planck equation oscillates rapidly with respect to 
the phase (j) = <f>\ + 4>2- We can average the probability 
density over this rapid phase |9j in the leading approxi- 
mation over 1/kL. This averaging can be most efficiently 
done on the level of the Fokker-Planck equation and, as 
the result, we obtain a partial differential equation on the 
density P(j], T\,T2) where i] — <pi — cf>2 is a "slow" vari- 
able. This is done, it is more convenient to come back 
to the level of the Langevin-type equations for only three 
variables 77, t\, and t%. Such a procedure requires in- 
troducing of two independent random forces fi(x), f2(x) 
rather then one. The reduced set of the Langevin type 
equations is found to be 



V 

Tl 
T 2 



1 sin 77 cos 77 + sin rjfi (x) 

1 



\ +ti (cos 77/1(2;) -sin77/ 2 (x) - ± cos 2 77) (10) 

\ . p I \ T 9 \ 



\ +T 2 {cosr)f 1 (x) 



- sin 77/2 (x) — I cos- 



n) 



The first equation is of the most importance and gov- 
erns the evolution of the slow phase 77. It contains the 
only random force /1 : (/i(x)/i(y)) = (2/£)5(x — y). The 
other two equations depend on additional random force 
/ 2 , which is independent of fi and has the same variance 
2/£. In the initial point x = the variables 77, n, and r 2 
have to be equal zero. The value of the parameter q in- 
fluences the dynamics of the phase 77: when q is positive 
(negative) the function rj(x) is also positive (negative). 
The limit q = leads to the trivial stationary solution 
r)(x) = and two equations for t\, r 2 become equivalent 
and can be used for calculating the time delay distribu- 
tion ||. In a general case q 7^ the last two equations 
still can be solved: 



1 



dx 



,f x dx'(fi cosrjT/2 sinr;-i cos 2 r/) 
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To perform the disorder averaging we have to integrate 
over the random forces /i, fi with the Gaussian measure. 
When q = the averaging of the product tit 2 is readily 
done and leads to the second moment of the time delay 
distribution: 



1 



2A 2 



(e 2T - 2T - 1) ; 



T = — 



(12) 



which increases exponentially with the length of the sam- 
ple. It is even simpler to check that (r) = L/k ~ A^T. 
When q ^ the averaging of the product tit 2 over the 
random forces gives rise to the correlation function (^) 
where u — 2q^. The integration over the noise f% is 
still Gaussian and can be easily done. To perform the 
other integration we restrict our consideration to the in- 
terval r\ £ (0, 7r) because of the obvious periodicity of 
the solutions and take advantage of the new variable 
g = In (tan 7^/2) to rewrite the first of Eqs.(^0|) in the 
following form: 

fi=9-V(g); V(g) = 2q cosh g + C'tanhg (13) 

This gives a possibility to replace the integration over fi 
by that over g taking into account the initial condition 
g(0) = —00 and the corresponding Jacobian: 



J = exp 




(14) 



The next step is to make use of the formal analogy with 
the Feynmann path integral of the quantum mechanics. 
After straightforward calculations we come back to the 
variable 77 and arrive at the time delay correlation func- 
tion (0) in the following form: 



K{u, T) = 2 f dh f 1 dt 2 R{T, h,t 2 ;uj)- T 2 
Jo Jo 



R = Cdr) S^y^e^-^ce^-^"" 1 



H L = sin 2 77 d 



H c = 
Hr = 



drj 

dry ^ 'H drj 

^2 . o 



/ 1 d 



(15) 



(16) 



'dri> 



sin rj 



ui — 

drj 



The action of the differential operators exp Hl.c,r a P~ 
plied to the right-hand side is equivalent to finding the 
solution of the corresponding evolution equations. Un- 
fortunately the eigenstates and eigenfunctions of the op- 
erators Hl,c,r are poorly understood and an exact cal- 
culation of the quantity R for arbitrary sample length 
L and frequency u> is beyond our possibilities. Still, the 
approximate solution can be found. It is worth notic- 
ing that the equations of this kind are common in the 
context of one dimensional localization and have been 
disscussed in the literature The thor- 

ough analysis of Eqs.(pl||l6|) shows that for exponentially 



small frequencies u> <C exp (— T) the main contribution to 
the integral over t\, ti in Eq.(|l5|) comes from the region 
t\ ~ <2 "C T. In this situation the time delay correla- 
tion function is proportional to lim J) _,o exp (THr)1. The 
latter expression is appropriate for the perturbation anal- 
ysis in u> which allows to present the result in the form 
of the asymptotic series: 



k=0 



2fc e (2fe+l)(2fc+2)T 



(17) 



where 7fc does not depend on T. With the help of the 
relation (^) we obtain j6|,[7|JTl|] that the resonance widths 
are distributed according to the log-normal law in the 
region of exponentially narrow resonances y <C exp(— T): 



p(y) = N^y- 3 / 2 exp (-(AT)- 1 In 2 y) 



(18) 



in agreement with natural expectations [|3|j7|,|l IB ■ 

Eqs.( |l5| , ^6|) can be analyzed also in the thermodynamic 
limit T — > 00 or large enough frequencies exp(— T) <C oj. 
The main contribution to the integral Eq.(|15|) exactly 
cancels the term T 2 and the leading correction turns out 
to be independent on T. 



K{u) = 



2 8 



d M 2 r(^) (|^|/8)^ 
M r 2 (1 + f ) cosh 2 ?f 



(19) 



The result of the numerical evaluation of the present in- 
tegral is shown in Fig.l. 

In the derivation of the expression above we con- 
sider the operators Hr ■ c,h acting in the Fourier space 
exp(2im?7). The drastic simplifications occur when we 
let the index m to be continous which is the correct ap- 
proximation for \m\ 3> 1. Such a limit should be ap- 
propriately matched with that for small m: \mu\ <C 1 
]lT| p"3| . The continous approximation is responsible for 
an unphysical tail for large frequencies. We expect, how- 
ever, that Eq.([i9|) is valid for u < 1. It is interesting to 
note that even though a simple one-pole approximation 
of the integral (|l9| ) gives rise to u^ 1 dependence in ui — > 
limit, such a dependence takes place only for extremely 
small frequencies: u> < 10~ 8 . The best fit for the correla- 
tion function in a very broad frequency region is oj -125 . 
Let us derive now the resonance width distribution p(y) 
which follows from Eqs.fl|l9|). Taking into account the 
identity 



dy y a+1 {y 2 +U 2 )- 1 = -7r|w| Q /(2sin7ra/2) (20) 



we restore the resonance density in the following form: 



p(y) 



n 



dp 



p 2 T 



' 3-H/J x 



sin 



7r(ifi— 3) 



r 2 (1 + f ) cosh 2 ffi 



(21) 
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This integral shows basically the same features as 
Eq. (|l9|) , (see Fig.||). The resonance widths turn out to 
be virtually cut at T ~ A £ /8 and the rest part of the plot 
should not to be taken seriously being an artifact of the 
approximation used. 
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FIG. 1. The time delay correlation function as given by 
Eq.(|l^) (solid line). The negative part of the function shown 
in the insert is the artifact of the approximation done. The 
dotted line corresponds to the dependence lu~ 1,25 . The func- 
tion dropes down to zero at the point to ~ 1/8. 



can argue that an effect of the open edge can be effec- 
tively simulated by adding the imaginary shift 7 = in 
to last diagonal element of the matrix. We picked up 
eigenvalues from the center of the spectrum and inves- 
tigated statistics of their imaginary parts. As shown in 
the Fig.|| the numerical results agree reasonably with our 
analytical predictions. 

In conclusion we address analytically the time de- 
lay correlations and resonances for the problem of re- 
flection from ID disordered sample. We find the reso- 
nance density p(T) to reveal the log-normal behaviour 
for the exponentially small widths and the algebraic de- 
pendence close to r~ 1,25 in the wide parametric range 
exp(— L/£) <C T/Aj <C 1. The time delay correlations 
are found to demonstrate similar behavior as a function 
of frequency. 
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FIG. 2. The resonance widths distribution calculated by 
the use of Eq.(|l|) (solid line) as compared with the results 
of the numerical simulations (squares) done for the matrices 
of the size 300 * 300. The insert and the dotted line are like 
those in Fig.[j]. 

To check our results we considered the simplest random 
matrix model which was expected to belong to the same 
universality class. Namely we diagonalize numerically as 
many as 10000 tridiagonal matricies of the size 300 * 300 
with its diagonal elements being uniformly distributed in 
the interval (—1,1). The off-diagonal elements are cho- 
sen to be equal to 1. In the same manner as in |l| we 
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